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ABSTRACT 


This  report  presents  the  solution  of  a  problem  involving  the  towing  of 
a  hydrofoil -cable  system  in  a  non-uniform  stream.  The  system  consists 
of  a  weightless  hydrofoil  connected  to  the  towing  vehicle  by  a  weightless 
flexible  cable.  The  towing  vehicle  moves  in  a  straight  line  at  constant 
speed  through  a  fluid  medium  wherein  the  steady  fluid  velocity  is  parallel 
to  the  direction  of  motion  but  varies  in  magnitude  in  the  normal  direction. 
The  problem  consists  in  determining  the  system  configuration  and  the 
cable  tension  for  specified  velocity  conditions,  cable  dimensions  and 
hydrofoil  characteristics.  It  is  shown  that  the  solution  of  this  two- 
dimensional  problem  involves  the  simultaneous  solution  of  two  trans¬ 
cendental  equations,  thus  requiring  the  use  of  numerical  and/or  graphical 
methods.  Finally,  a  numerical  example  is  presented  to  illustrate  the 
method. 
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NOTATION 

A.  PHYSICAL  QUANTITIES 

A  Representative  (planform)  area  of  hydrofoil 
a  Ratio  of  stream  vel<  city  gradient  to  velocity  of  towing  vehicle 

b  Fluid  velocity  gradient  (constant  in  numerical  example) 

c  Constant  (  =  k/T' ) 

Cj^  Hydrofoil  drag  coefficient 

CL  Hydrofoil  lift  coefficient 

CR  Cable  resistance  coefficient  (for  cable  normal  to  stream) 

D  Hydrofoil  drag  force 

P  Cable  diameter 

c 

2 

k  Reference  normal  cable  loading  (  =  CR  •  4-  V*  Dc> 

L  Hydrofoil  lift  force 

N  Cable  normal  force 

S  Total  cable  length 

s  Running  cable  length 

T  Cable  tension 

T'  Reference  tension  value  ( =  y/cL  2  +  CD  2  .  A  •  Vo  2  ) 

VQ  Velocity  of  towing  vehicle,  relative  to  still  water 

V8  Velocity  of  stream,  relative  to  still  water 

x,  y  Rectangular  coordinates  of  cable  element 
P  Mass  density  of  fluid 

4>  Local  inclination  of  cable  curve 
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B.  MATHEMATICAL  QUANTITIES 

See  Eq.  (15) 

See  Eq.  (7) 

See  Eq.  (25) 

See  Eq.  (27) 

See  Eq.  (22) 

e  -cF<yi)/f(yi) 

C.  SUBSCRIPTS 

o  Junction  of  cable  and  towing  vehicle  (origin  of  coordinates) 

1  Junction  of  cable  and  hydrofoil 

a  Denotes  value  corresponding  to  constant  value  of 

b  Denotes  value  corresponding  to  constant  value  of 
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SECTION  1 
INTRODUCTION 


1.  1  BACKGROUND 

Under  ONR  Contract  N0nr-3582(00),  a  general  analyeie  is  to  be  made 
of  the  steady  turning  of  cable -towed  body  systems,  wherein  the  towing 
vehicle  moves  in  a  circular  path  at  constant  speed. 

This  is  a  relatively  complex  problem,  involving  several  features  never 
previously  analyzed,  as  follows: 

(1)  A  non-uniform  field  of  flow; 

(2)  Motions  and  forces  in  three  dimensions; 

(3)  The  complicated  relation  between  the  motion  of  the  towed  body  (as 
part  of  the  system),  its  attitude,  and  the  forces  and  moments  acting  upon  it. 

Accordingly,  it  was  considered  desirable  to  analyze  a  number  of  non¬ 
trivial  problems  involving  one  or  more  of  these  features  in  order  to 
establish  mathematical  procedures  for,  and  gain  insight  into,  the  physical 
aspects  of  the  general  turning  problem. 

This  report  covers  an  analysis  of  the  first  such  problem  established, 
involving  a  non-uniform  flow  field.  In  addition  to  its  pertinence  to  the 
turning  problem,  the  particular  case  analyzed  herein  also  has  appreciable 
intrinsic  interest. 

1.  2  STATEMENT  OF  PROBLEM 

A  neutrally  buoyant,  completely  submerged,  smooth,  round  cable  of 
specified  diameter  and  length  is  used  to  tow  a  stable  hydrofoil,  as  illustrated 
in  figure  1.  The  hydrofoil  and  its  attachment  to  the  tow  cable  are  arranged 
so  that  the  hydrofoil  operates  at  a  fixed  angle  of  attack  and,  consequently, 
at  fixed  (specified)  lift  and  drag  coefficients.  The  cable  is  towed  by  a 
vehicle  which  moves  horizontally  in  a  straight  line  (x-axis,  negative  x- 
direction)  at  a  specified  constant  speed  through  a  moving  fluid.  The  fluid 
velocity  is  everywhere  horizontal  and  its  magnitude  is  a  specified  function 
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of  the  y-coordinate.  The  problem  then  is  to  establish  a  procedure  for 
determining  the  hydrofoil  location  relative  to  the  towing  vehicle,  and  the 
cable  shape  and  (constant)  tension,  if  the  cable  is  subject  only  to  a  normal 
force  per  unit  length  proportional  to  the  square  of  the  local  normal  component 
of  relative  velocity  between  cable  and  fluid. 
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SECTION  2 

METHOD  OF  SOLUTION 

2.  1  ANALYSIS 

For  system  equilibrium,  the  entire  cable  and  the  hydrofoil  must  move 
with  constant  velocity,  Vo,  (relative  to  an  inertial  frame,  or  "still  water") 
in  the  same  direction  as  the  towing  vehicle.  As  shown  in  figure  1,  this 
direction  is  taken  as  the  negative  x-direction.  The  (specified)  stream 
velocity  (relative  to  still  water),  taken  as  positive  in  the  x-direction,  can 
be  written  as  Vg(y).  Hence,  the  fluid  velocity  relative  to  the  hydrofoil- 
cable  system  is: 

V  (y)  =  Vq  +  Vs  (y)  (1) 

It  is  assumed  first  that  the  only  hydrodynamic  forces  acting  on  the 
cable  are  normal  forces,  i.  e.  the  frictional  (tangential)  forces  are  neglected. 
Further,  the  single  equation  defining  the  equilibrium  of  the  cable  element  is 
given  by  Eq.  (2)  of  Reference  1.  In  the  present  notation,  this  is: 

T  d<t>  =  dN  (2) 

where:  T  =  Cable  tension  (unknown  constant) 

4>  -  Angle  of  inclination  of  cable  element  (to  x-axis) 

dN.  =  Normal  force  on  cable  element 

In  accordance  with  Reference  1,  the  hydrodynamic  force,  dN,  is  taken 
as  proportional  to  the  square  of  the  normal  component  of  the  local  relative 
fluid  velocity,  specifically: 

dN  =  CR  •  (V  sin  *  ) 2  Dc  ds  (3) 
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Combining  Eq.  (2)  and  (3): 


dN 

d<£ 


V2  »in 2  4> 


ds  _ 


“17 


=  T 


(4) 


which  can  be  written  as: 

k  •  f(y)  •  sin  2  4>  •  =  T 


where: 


k 


CR- 


(5) 

(6) 


and 


v«  (y> 


(?) 


i.  e. ,  f  (y)  is  a  known  function. 

It  is  seen  that  Eq.  (5)  involves  two  (2) unknown  functions, y  (0)  and  s  {<P)t 
as  well  as  the  unknown  constant,  T. 

Now,  the  tension  must  equal  the  resultant  hydrofoil  load: 

T  =  /l2  +  D2 

-  W=L2  ♦  cD2  f'-f-  v,2>) 

wherein  A,  CL  and  are  known  quantities,  but  Vj  =  V  (yj)  is  unknown. 
Eq.  (8)  can  be  written  as: 


T  =  T'  .  f(yj) 


(9) 
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where  T  1  is  the  known  constant: 

T '  =  A 


/  cT  2  + 


L  T  CD 


r 

+ 


(10) 


and  is  the  unknown  total  cable  height. 

Combining  Eq.  (5)  and  (9),  and  using  the  relation: 
dy  =  (  sin  ^  )  ds, 


there  is  obtained: 


d  4> 

sin  0 


=  c  • 


Mr  * 


where  c  is  the  known  constant: 


c  = 


1 


CL  +  CD 


Taking  the  origin  of  coordinates  (and  of  cable  length)  at  the  towing 
vehicle,  as  shown  in  figure  1,  integration  of  Eq.  (12)  yields: 


In 


4>  , 

~T  \ 


tan  -2-  > 


tan 


JLM. 


J 


*  (Vi) 


(H) 


(12) 


(13) 


(14) 


where  4>Q  is  the  cable  inclination  at  the  towing  vehicle,  and  F  (y)  is  the 
integral: 

F(y)  *  *(y)  dy  (15) 

As  f(y)  is  a  known  function,  so  is  F(y). 

At  the  hydrofoil -cable  junction: 


4>  =  <t> 


i  •  *  9  *i 


(16) 
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Hence, 


-Jl  c‘ 
■+[■ 


esc  ^  dy 


(  tan  -y-2- )  Gj  (y,  yj )  +  (cot  — ~ 


■i  <v>  vi>] 


where: 


0,(y.  y,)  ./oy  g  dy  -/J  .cr<y>/f(y,>  dy 
°-i  IV.  v,)  =/: .  •>  dy ./;  .-cF(v>"<yi>  dy 


Finally,  the  total  cable  length  is  given  by: 


=  -j-  (tan  ).  T j  (yj)  +  (cot  ) 


where  (y  j  )  and  T  ^  (y^  )  are  the  known  functions: 

ri  Ivi>  *  o,  (y, .  y,)  -J*»  e  cF<y)/,(yi’  »y 


r..  <y.)  *  °.,<y,.yi)  *J?1  •  *cF(y,/f<y>)  dy  ) 


Consequently,  for  a  given  value  of  S,  Eq,  (26)  represents  another 
definite  functional,  or  graphical,  relation  between  y^  and  4 Q : 

V  *o,b'Vl> 
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The  simultaneous  solution  of  Eq.  (20)  and  (28),  which  must  generally 
be  obtained  by  numerical  and/or  graphical  methods,  furnishes  the  solution 
of  the  present  problem,  i.  e.  permits  determination  of  ^  Q  and  y  j  for 
prescribed  values  of  ^  ^  and  S.  The  cable  tension  can  then  be  found  from 
Eq.  (11). 

In  many  cases,  it  is  also  desirable  to  determine  the  shape  of  the  cable. 
This  curve  is,  of  course,  given  in  intrinsic  form,  s  =  s  (^),  by  Eq.  (21) 
and  (24).  Its  Cartesian  equation  is  readily  shown  to  be: 

,  r  *  *  "i 

X  =  -L-  jjcot  )  G_J  (y,  yt)  -  (tan -^2. )  Gj  (y,  yt)j  (29) 

Similarly,  the  horizontal  distance  between  the  hydrofoil  and  the  towing 
vehicle  is : 

*,  =  -j-  [(cot -la )  r_,  (y,)  .(«»n^)  T,  (Vl )  J  (JO) 
2.2  COMPUTATION  PROCEDURE 

Numerical  calculations  can  be  performed  using  the  following  procedure, 
wherein  the  choice  of  analytical,  numerical  and/or  graphic  methods  for 
each  step  is  left  to  the  discretion  of  the  reader: 

A.  Using  Eq.  (10),  (13),  and  (18),  calculate 
T  '  ,  c,  and  4>  j. 

B.  Combine  the  towing  velocity,  Vo,  with  the  known  fluid 
velocity  profile,  V  (y),  (with  proper  signs),  to  obtain: 


which  is  functionally  equivalent  to  f  (y^  ). 
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C.  Calculate  the  function: 

F  (y)  =  J**  i  (y)  dy, 

which  is  functionally  equivalent  to  F  (y^  ). 

D.  Calculate  the  function: 

<yi>  * 


E.  Using  the  value  of  ^  ^  from  Step  A,  plot  the  curve  of  ^  Q  ve.  yj 

corresponding  to  Eq.  (19): 

4>q  =  2  tan  4  |\tan  (Yj  >1 

F.  Calculate  the  two  families  of  functions: 


g  (y.  = 

and  its  reciprocal. 


e  c  F  (y)  /  f  (yj) 

g  _1  (y.  y^ 


G.  Calculate  the  two  functions: 

ri  (yi 1  *  JJ1  «  (y-  yi  >  dy 

r-i  <yi>  *  (y>  yi>  dy 


H.  Using  the  known  value  of  S,  use  Eq.  (26)  to  plot  vs.  y 
that  for  fixed  ,  Eq.  (26)  is  a  quadratic  equation  for  tan 


I* 


(Note 

). 


* 
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I.  Superpose  the  two  curvet  of  vs.  obtained  in  Steps  E  and  H. 
Their  point  of  intersection  determines  the  desired  values  of  and 

V 

J.  Using  the  latter  value  of  y^  ,  determine  the  numerical  value  of 

f  (y^ )  from  Eq.  (7)  and  thus  the  cable  tension.  T,  from  Eq.  (9). 

K.  To  find  x  ^  ,  use  the  final  values  of  y  ^  and  4> Q  in  Eq.  (30). 

L.  To  find  the  cable  curve,  use  the  final  value  of  y^  to  calculate  the 
functions  G ^  and  G  ^  defined  by  Eq.  (25)  and  combine  these 
according  to  Eq.  (29),  using  the  final  value  of  ^  . 


2.3  NUMERICAL  EXAMPLE 

To  illustrate  the  preceding  analysis,  calculations  were  made  for  a 
hydrofoil -cable  system  towed  through  a  fluid  having  a  constant  velocity 
gradient,  i.  e.  a  linear  velocity  distribution.  Three  different  gradient 
values  were  assumed  while  all  other  system  characteristics,  including 
towing  speed,  were  kept  constant,  as  follows: 


Cable: 
Hydrofoil: 
Medium: 
Towing  Speed: 
Fluid  Velocity: 


CR  =  1.  2,  D  =  1  in. ,  S  =  400  ft. 

CL  A  =  20  ft.  2,  CD  A  =  3-1/3  ft.2 

Sea  Water,  P  =  2  slugs /ft.  * 

V  =5  knots 
o 

Vg  =  by,  b  =  0,  0.5,  1.  0  knots/ 100  ft. 


The  cable  and  hydrofoil  characteristics  yield  the  values: 

c  =  0.493  per  100  ft. 

<t>l  =  80*  32' 
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and  the  basic  functions  are: 

f(y)  =  (1  +  ay)2, 

F  (y)  =  y  (  1  +  ay 

Following  the  procedure  outlined  above,  and  using  numerical  methods 
in  conjunction  with  an  IBM  610  digital  computer  to  evaluate  all  of  the 
required  integrals,  the  pair  of  curves  (  4>  Q  vs.  y ^  )  for  each  gradient 
value  was  obtained.  All  of  these  are  shown  in  figure  2. 

The  resulting  solutions  for  4>  and  y^  ,  as  well  as  the  corresponding 
values  of  and  T  have  been  plotted  against  the  gradient  parameter,  b, 
in  figure  3.  It  is  seen  that  for  the  assumed  system,  the  effect  of  this  type 
of  velocity  distribution  is  to  elevate  the  hydrofoil  with  respect  to  the  towing 
vehicle,  thus  increasing  the  cable  tension.  The  shapes  of  the  cable  curves 
were  not  calculated. 


a  =  b/V  =  0,  0.  1,  0.  2  per  100  ft. 
o 
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SECTION  3 
DISCUSSION 

The  following  facts  may  now  be  noted: 

1.  UNIFORM  STREAM 

It  is  readily  shown  that  the  present  method  correctly  yields  the 
known  solution  for  the  case  of  the  uniform  stream,  wherein  the  cable  has 
the  shape  of  a  catenary  curve.  (Reference  1.) 

2.  RELATED  PROBLEMS 

By  minor  variations  of  the  present  method,  i.  e.  by  use  of  different 
types  of  cross -plots,  the  relations  established  above  can  be  used  to  solve 
any  problem  wherein  any  two  of  the  four  quantities:  4>Q  ,  Q  j  ,  y»  ,  S  , 
have  to  be  determined  from  known  values  of  the  other  two. 

3.  MORE  GENERAL  CASES 

The  present  (i.  e.  cross -plotting)  method  can  be  extended  readily  to 
solve  more  general  cases  such  as  those  involving  fluid  friction  and  finite 
cable  weight.  Evidently,  however,  the  labor  involved  in  obtaining  numerical 
solutions  increases  rapidly  as  the  physical  system  becomes  more  complex. 
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SECTION  4 
CONCLUSIONS 


1.  For  the  two-dimensional  hydrofoil -cable  system,  towed  in  a 
straight  line  at  constant  speed  in  a  non-uniform  stream,  the  towing  problem 
cannot  generally  be  solved  in  an  explicit  manner,  even  when  only  normal 
hydrodynamic  forces  on  the  cable  and  weightless  elements  are  assumed. 
Instead,  the  problem  has  to  be  solved  by  an  indirect  method,  i.  e. ,  by 
simultaneous  (numerical,  graphical,  etc. )  solution  of  two  transcendental 
equations. 

2.  The  labor  involved  in  obtaining  numerical  solutions  is  appreciable 
and  would  be  substantially  increased  in  the  more  general  cases  involving 
fluid  friction  and  gravity  forces. 

3.  With  regard  to  the  ultimate  problem  involving  steady  turning  in 
three  dimensions,  it  is  concluded  that  those  aspects  relating  to  the  non- 
uniform  flow  condition  do  not  offer  any  inherent  mathematical  difficulty, 

but  will  certainly  add  considerable  labor  in  obtaining  numerical  solutions. 
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(o)  VELOCITIES  RELATIVE  TO  STILL  WATER 


(b)  VELOCITIES  RELATIVE  TO  HYDROFOIL-  CABLE  SYSTEM 


FIGURE  I.  SYSTEM  CONFIGURATION,  ILLUSTRATING 

ABSOLUTE  AND  RELATIVE  VELOCITY  RELATIONS 


FIGURE  3.  SYSTEM  GEOMETRY  AND  TENSION  VALUES 

V  S 

GRADIENT  MAGNITUDE 


